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Abstract 

We point out that adding higher dimension operators to the Nambu- Jona- 
Lasinio model can stabihze chiral symmetry breaking and resolve the issue raised 
by Kleinert and Van den Bossche that chiral symmetry breaking cannot occur for 

^ ^c"^ due to large quantum fluctuations of the composite fields. 
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The Nambu-Jona-Lasinio (NJL) model has been an important tool in studying chiral 
symmetry breaking [|T], ^, ^. The model is simple, comprised of four-fermion contact 
interactions, and solvable in the large Nc (number of colors) limit, showing chiral sym- 
metry breaking in a very clean manner [Q. Recently, however, Kleinert and Van den 
Bossche (KV) argued that the strong quantum fluctuations of the mesons around the 
(classically) broken vacuum in the low energy effective Lagrangian of the NJL model could 
restore chiral symmetry. This chiral symmetry restoration occurs when the so-called stiff- 
ness, which is the pion decay constant squared in the NJL model, becomes smaller than 
a certain critical value. This happens in the real world case Nc = ?> and can cause a 
serious problem for using the NJL model as a low energy effective description of the chiral 
symmetry breaking in QCD. 

The four-fermion contact interactions of the NJL model, which makes the model non- 
renormalizable, may be thought as effective interactions arising from a more fundamental 
one, for example from an exchange of massive intermediate particles in a renormalizable 
theory J^, In this case, generally there would be higher dimension operators (HDO's) 
in addition to the contact interactions. These operators will be suppressed by the cutoff 
of the model, and at tree level have small effects and can be ignored. However, they 
cannot be ignored at loop level as long as momentum cutoff is used in regularization. 
This becomes clear when we recall that some of the loop diagrams in the NJL model 
are quadratically divergent, and thus in momentum cutoff the divergence could cancel 
the cutoff suppression of the HDO's, giving unsuppressed contribution to the low energy 
effective Lagrangian. 

The importance of HDO's in the NJL model was previously emphasized by several 
authors [|, 0, while it was first discussed by Suzuki in a specific form of dim=8 



operators and by Hasenfratz et al. [|T2| in a more generalized form in the context of top 
quark condensation model of electroweak symmetry breaking. In this note we investigate 
the effects of HDO's on the chiral symmetry breaking and point out in particular that 
adding proper HDO's to the NJL model could suppress the quantum fluctuations of the 
mesons and prevent chiral symmetry restoration. 
We start with a two flavor NJL model: 



where the massless quarks = {u, d)\ i = 1, ■ • ■ , carry Nc colors, are the Pauli 
matrices, g is a dimensionless coupling and A is the cutoff. Introducing auxiliary scalar 
(cr) and pseudoscalar (tTq) fields the Lagrangian can be written as 

NcA"^ 

C = - a - i'j5TaTTa)i' ^ (cr^ + TT^). (2) 

Integrating out the quark fields we obtain an 0(4) symmetric effective Lagrangian for the 
mesons 



where = + vr^ , while each term in C has at least four derivatives and inverse power 
suppressed in the cutoff A^, thus not of our interest. The wave function normalization 
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Figure 1: Wave function normalization of the mesons by a quark-antiquark loop. 



Z{p^) arises from Fig. 1 and is given by 
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while the potential V reads 
V^p") = 



P 



(5) 



where pe, Qe denote the Wick rotated Euclidean momenta. Throughout the article we 
compTitc as usual in the large Nc limit. When the coupling g is sufficiently large the 
potential can have a minimum at a nonvanishing = and the symmetry is broken in 
the large limit. The minimum of the potential is determined by the gap equation 
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which has a nonvanishing solution provided g > Qcr, where g^ — 27r^. 

Expanding about the vacuum (a, Tr^) = (po, 0), and putting a = a' + po, we obtain the 
quadratic part of the Lagrangian (the prime in a' is now omitted) : 



Z{pl) 



{d^ay - Apy + {d^T^a) 



+ 



(7) 



which shows that cr has a mass = 2po and tTq are the massless Nambu-Goldstone 
bosons. 
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KV observed that although the potential (|^) breaks the chiral symmetry in the large 
Nc limit, the quantum fluctuations of the mesons can restore the symmetry when Nc is 
smaller than a certain critical value N^"^. This observation was made as follows. The 
nonlinear sigma model of the massless modes associated with the Lagrangian is given 
by 



c 



where S = (cr, TCa) and A(x) is the Lagrange multiplier for the constraint S S — pg = 
Integrating out the T,^ gives an effective function in A 



Wi\) 



(9) 



where /? = Z{pl)pl is the stiffness. Extremizing the function a second gap equation for A 
is obtained: 

— ^ / dt 

87r2 [Jo t + X 



(10) 



which has a nonvanishing solution A = Aq when the stiffness satisfies 



P<Pc 



(11) 



where Pc 
mass m" 



A'^/Att'^. a nonvanishing solution for A implies that the obtain a common 
= Aq and, consequently, the chiral symmetry remains unbroken. Using @ it 
can be checked that for A^^ = 3 the condition is satisfied when p^/A^ < e~^, which 
implies no chiral symmetry breaking in the real world case. Here an important note is in 
order. In drawing the conclusion it was assumed that the UV cutoff in the mesonic loop 
(in eq.(P)), which is 1/Nc correction, was identical to that of the quark loops. Since the 
theory is nonrenormalizable there is, however, no reason that the assumption should be 
true, and in fact the cutoffs should be assumed different |]l3l . A recent study shows there 
is a possibility for the chiral symmetry restoration only when the cutoff for meson loops 
is not small compared to that of the quark loops p^ . 

From the discussion so far it is clear that the problem lies with too small a wave 
function normalization Z{pq) for a given po- We now suggest a modification of the NJL 
model in which the stiffness can be sufficiently large so that (3 > /3cr even for = 3. 
Back to (0), we see the easiest way to increase the wave function normalization is to add 



cNr 



(12) 



to the original Lagrangian, where c is a positive constant. This new term obviously shifts 
the wave function normalization by the amount 



6Z = cNr 



(13) 



Since the addition of (|^) does not affect the potential V{p ) the gap equation is not 



modified, and the shift ( 1131) can render the stiffness greater than the critical value when 
c is sufficiently large. Numerically one can obtain N^"^ < 3 with c = 0.7 /g. 
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The addition of the term (^2[) corresponds to an addition of HDO's to the original 
NJL Lagrangian ([^). To see this we may solve iteratively the equations of motion for the 
meson fields arising from the Lagrangian (0) and (p!2D, obtaining 



a 



(1 - cgd^/K^ + • ■ OV-V^, 



+ cgd^ / ^'r'i'l,T,i' 
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which leads to a modified NJL Lagrangian 
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(1 - cgd^/K^ + • ■ OV-Tsra^, 
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This clearly shows that the wave function normalization is sensitive to the presence of 
HDO's. 

We may now consider a more general modification of the NJL model: 
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where 



K{d'^IK\ 9VA') = 1 + (s'Va^) (d^/x^y 

is a finite polynomial so that the Lagrangian is local, while cio, coi, 
constants. Here 9'^ =d^d^ and 



(17) 

are dimensionless 
(18) 



where Cj are dimensionless constants and the parameter m is assumed to be rn^ ^ p^. The 
derivatives in K are defined to act on the quark fields only. 

Integration over the quark fields gives an effective Lagrangian of the mesons in the 
form (H), now with the wave function normalization and the potential given by 
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Minimization of the potential gives the gap equation as 
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Using this equation and 
the HDO's as 



where 



we can now obtain the shift in the critical coupling due to 
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(24) 



which shows that the critical coupling depends on the HDO's, unsuppressed in l/A^, as 
it should be since the integral in the gap equation is quadratically divergent. 
Now the wave function normalization in (p^Of ) can be written as 

(25) 



Zip') = Zoip') + /i.(Q„ci) + o(pVa'), 



where Zo(p^), the wave function normalization of the original NJL model, is given by (^, 
and hz is a function depending only on the coefficients of the HDO's. This can be easily 
shown by noticing that the coefficient of the logarithmically divergent term ln(A2/p2) in 
Z(p2) does not depend on the presence of the HDO's, which can be seen by expanding 
the integrand in (|T9|) in powers of the coefficients Cij and performing the integration. This 
shows that the effect of HDO's in K to the wave function normalization is not suppressed 
by 1/A2 but 0(1). 

Expanding the potential (^) about the vacuum we obtain the a mass 



where 
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(27) 



As in the case of Z{p''), the coefficient of the logarithmically divergent term in the integral 
is independent of the HDO's, and thus (^) may be written as 



V'ipl) = Zoipl) + hm{c,,. C2) + 0{pl/A' 



(28) 
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where hm{cij, C2) is a function of Cy, C2 only. Using p5|) , (p^) we have 

= 2po^ Y+'^hjll [1 + ^(Po/A^)] • (29) 

This shows that with the presence of HDO's the standard relation rria = 2po is no longer 
valid [|ll|]. Only when the coupling g is fine tuned so that Po/A'^ <^ 1, is about twice 
the constituent quark mass. In fact this is just one evidence for a more general conclusion. 
The logarithmic divergence in (|l^), (|27|) and in the coefficient of the marginal operator 
from the expansion of the potential ( pT]) shows that the effect of HDO's disappears 
logarithmically in p^/A^ from the low energy effective Lagrangian. Of course this conclu- 
sion is derived from the large Nc calculation, and when Po/A^ <C 1 the renormalization 
group running of the parameters in the low energy effective Lagrangian due to the meson 
fluctuations should be considered as in the top quark condensation model |TB|, |T^. 



In the NJL model for QCD, however, Pg/A^ is not so small (po ~ 350 MeV and A ~ 900 
MeV), and consequently HDO's can play a role in determining the phase of the system. 

In conclusion we have introduced HDO's in the NJL model through the auxiliary 
meson fields and shown that they can have significant effects on the low energy effective 
Lagrangian when the NJL coupling is such that p^/A"^ <^1. In such case the low energy 
quantities such as the stiffness or sigma meson mass are sensitive to the presence of HDO's 
and should be regarded as free parameters unless sufficient informations are available to 
determine the form and magnitude of the HDO's. This then allows one to avoid the chiral 
symmetry restoration suggested by KV by modifying the NJL model with proper HDO's. 
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